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Abstract 

The scale factors of an arbitrary orthogonal space are a measure of its content of 
homogeneous orthogonal space. In the present study, it is shown, that their spatial and 
temporal rates of variation do not contribute to the differential calculus of arbitrary 
functions in orthogonal space. Based on this, the Navier-Stokes equations are formu- 
lated accordingly, to provide a method of studying the kinematics of fluid motion in 
orthogonal plane space. Employing this formulation and regular perturbation theory, 
the kinematic physical measures of the flow of an incompressible, viscous fluid around 
two identical circular cylinders, which move with equal and opposite velocity towards 
central impact are evaluated. For this case, the space is parametrised according to the 
bipolar transformation of the cartesian plane. 

1 Introduction 

The expressions of the basic differen t ial op erators, according to an arbitrary orthogonal 



vector base, is known, see e.g. ISpain I ([1967). The scale factors, which characterise such 



a base, are, in general, functions of space and time. Their spatial rate of variation has 
been considered to contribute to the spatial rate of variation of arbitrary functions of the 
orthogonal space. The present study shows that this is not true. The significance of the scale 
factors is that they are a measure of the deformation of the space, which they characterise. 

Based on this result, which affects the differential calculus of functions in orthogonal 
space, the Navier-Stokes equations are formulated in an arbitrarily parametrised, orthogonal, 
plane space. An application of this formulation follows, in order to obtain an asymptotic 
series solution for the kinematic physical measures of the flow of a fluid around two identical 
circul ar cylinders, whch moy e with equal and opposite velocity, towards central impact. 

Bampalas fc Graham ( 2008f) studied the flow, based on a numerical solution of the 



Navier-Stokes equations in the plane. However, as the boundaries of the cylinders approach, 
mesh resolution effects affect the accuracy of the numerical solution. 

The next section presents the expressions of the basic differential operators in orthog- 
onal space, on which the subsequent analysis is based on. In sectionfSl the formulation of the 
Navier-Stokes equations in an arbitrarily parametrised, orthgonal, plane space is presented. 
In section^ the bipolar transformation of the cartesian plane is employed, to adapt the for- 
mulation to the case of the flow of a fluid around two circular cylinders and the asymptotic 
series solution is obtained. 
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2 The basic differential operators in an arbitrary or- 
thogonal space 

Consider the right-handed, orthogonal vector base Ij, i G N : i G [1,3], the spaces q, 
s G R 3 and time t G K : t > 0. Implying the summation convention, 

q:=ftlj , s(q;i) := Sj(q;t) Ij ; 

Si (q; i) are independent variables in s and qi are independent variables in q. 
The units of spatial measure of q and s are 



dq := {dqili G q : dqi — > 0} =/= , cfe := {cfejlj G s : cfej — > 0} 7^ . 

The variation of arbitrary, continuous functions of q or s along dq or cfe is assumed to be 
linear. 

Define a mapping 

i 1 - {s^q : dsjlj =s i (q;t)dt? i I i , Sj(q;t) G R : Si ^ 0} . 

(3(q; t) := Sjlj represents a density measure of s in q. For the special case q, s G R 2 , 5, = s. 
s is measured as 



/ du Ij = / 0i[q(u);f]dtt I< 
Jo Jo 



Consider an arbitrary, continuous scalar function f(s,t) G R, with continuous deriva- 
tives. The spatial r ate of variation o f / along Ij is described by the directional derivative 
/ Sj along Ij (see e.g. Krevszig (1999)), which is defined as 

, r / .n, .1 df r /[s(q;t) + dsi(q;i)Ii,t] -/[s(q;i),t] 
f Sl [s(q;t),t = —:= lim = 

dsi dsi^o dsi 

/[s(q; t) + Si(q; t)dq l l l ,t] - /[s(q; i), t] 



= lim 



s,(q;t)dg ! ^o Sj(q;t)dgj 

-l, ,s r /[q + jgJi, 6(q; Q] - /[q, t; 6(q; f)] 
= 5^ (q;i) hm 

d qi ->o c% 



Therefore, 



/ Si [s(q;t),t] = / Si [q,t;6(q;t)] =S- 1 (q;i)/ 9l [q,i;©(q;i)] • 

From the first fundamental theorem of calculus (see e.g. lApostoll ()l967f )h the indefinite 
integral of f Si along Ij, is 

/[s(q;t),t] = /[q,t;6(q;t)] = | ' /^[s^),^!* = | ' / g( [q( U ),i;6(q;t)](iu , 

where £ G s • Ij, ( G q • Ij. 

Space and time are considered to be independent sets. Therefore, the temporal rate 
of variation of /[s(q; t), t] is 
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ft s(q;t),t = := hm — = / t q, t; 6(q; <)J . 

at dt-vo at 

The indefinite integral of ft along time is 

/[s(q; t),i] = / / t [q,u;6(q;t)]<iu , 6 [0,t] . 

Starting from the fact that s, (q; i) are in dependent va riables in s and following the 
way usually employed in vector analysis, see e.g. ISpain I ([19670 . the expressions for the basic 
differential operators in s (see the appendix) are 



V:=d Si Ii , V- : l d s 
V 2 := <9«2 , V 2 := 'V 2 !,; 



Vx := [ m{t+2 ^d Sm 
V 4 := 3,4 + 



_m(i+l) q 

+ d e 2 



„ , V 4 := l V 4 I, ; 



(9 denotes partial differentiation, a left superscript of an operator indicates the vector com- 
ponent on which it operates and 



m(i) :— 



i h- >• i , i £ [1,3] 
2 1 y % 3 , i £ [4, 5] 



The differential operators in q are obtained by applying ./# on the corresponding 
operators in s. The operators in s,q £ E 2 are obtained by those in s,q £ K 3 by setting 



d S3 = 0, d q: 



0. 



3 A normalisation of the Navier-Stokes equations in or- 
thogonal plane space 

Consider the mechanics of fluid of density p and kinematic viscosity in plane space 
s £ M 2 and in time. Applying a mapping ^# : {s M> q}, such as that obtained by a 
conformal mapping of the cartesian plane in q, it is 



6(q;f) =S;(q;i)Ii =s(q;i)Ij . 

The fundamental physical measures, of fluid mechanics in the plane, are the stream- 
function \I> = ^13, the velocity U = t/jlj, i — 1,2, the vorticity uj = WI3 and the pressure 
P. The Navier-Stokes equations in the plane are 

D t V = -/) _1 VP + vV 2 \J , V-U = , or 

D t u} = vV 2 uj , V 2 * = -u) ; 

D t signifies the material derivative 

D t := d t + U • V . 

The kinematic relation between the velocity and the streamfunction for plane motion is 

U = V x * . 
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The physical measures are functions of the independent variables and the parameters as 



$[s(q;t),£;i/;p] =*[q,*;s(q;t);i/;p] , U[s(q; t), t; v; p] = t/;[q, t;s(q; t); v\ p] , 
u}[s(q;t),t;v;p] = uj[q,t;s{q;t);v; p] , P[s(q;t), t; v; p] = P[q, i;s(q;t); v\ p] . 

Consider the boundary dq of q to be rectangular and symmetrical with respect to 
q = 0. Set the lines qi = ±k, to confine q and define 1 = ijlj. Introduce a characteristic 
speed U s , a characteristic, cartesian length scale I and for s(q;t) with units of length, set 
s(q;i) to scale length in s and s(q;t)U~ 1 time. The characteristic length and time scales, 
are functions of space and time, but based on the analysis presented in section [2 they are 
parameters in the differential calculus of the physical measures of fluid motion in q. The 
scaled variables, parameters and measures are 



Qi ■= Qil^ 1 

Sj(q;i;s) := s;(q;i)s _1 
t:= tUsS^ 1 
* := ^(UsS)- 1 
P := P(pU^)- 1 



i(q; t) := si 1 
d§i := ksdqi 
Pe s := UsSv- 1 
U := UU- 1 
4> := §s 



Si(q;t) := Si(q;?)f _1 

dsi := lidqi 

Re t := Uslv' 1 

u := wsC/ s _1 
l 



ujs 



Introduce the scaled differential operators 



V :=sV 



, V-:=sV- 
V 2 :=s 2 V 2 , Df := d~ t + U • V 

and the scaled Navier-Stokes equations in the plane are 

= —VP + Pe^V^U , 
D~ t Cj = Pe^V 2 ^ 
U = V x * . 

The differential equation for 'J is 

P t -(V 2 *) = Pe^V 4 * 

and in explicit form 



Vx := sV x 
V 4 := s 4 V 4 



v-u = o 

V 2 * = -Q 



The functional dependence of the scaled physical measures on the scaled variables and the 
parameters is 



*(s,f;s;Pe 5 ) = *(q, t; s; Re s ; 1) 
L>(s, t; s; Re s ) = w(q, i; s; Pe 6 ; 1) 



, U(s, t;S;Re s ) = U(q, t; s; Pe 5 ; 1) , 
, P(s,f;s;Pe 5 ) = P(q, i;s; Pe 5 ; 1) . 
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4 The kinematics of the flow of a fluid around two cir- 
cular cylinders moving towards central impact 

Consider two identical, circular cylinders of radius I — R, immersed in incompressible fluid. 
The fluid is set into motion due to the motion of the cylinders, with constant and equal in 
magnitude, but opposite velocity U s towards central impact (see figure[T]). 2h(t) — 2ho—2U s t 
denotes the minimum distance between the boundaries of the cylinders, ho denotes half the 
initial di stance between th e cylin ders and [=pa(t), 0] are the foci of the bipolar transformation 
(see e.g. iMilne-Thomson I ((i960)), denote the unit vectors of the cartesian plane (x,y). 
The scaled parameters and space variables in figure Q] are 

-- 1 sinh- 1 {[e(2 + e)] 1 / 2 } 
y 2 +sinh 2 (^)]- 1 } 

,2/ /M-n 



e(t) = fiR' 1 




TO 


(*) ■= Mt) 


tt" 1 = n- 1 sinh" 1 (aR~ 1 ) 


T(x,y\To) := 


t{ttt' q ) 


1 = M 


_1 tanh -1 {2 sinh(7rrQ)x[5; 




an 


1 


= 7T _1 tan - 


1 {2sinh(^)y[i 2 +y 2 - 


i(^,T;TQ) := 


xR~ 


-1 


= sinh(7TTQ 


) sinh(7rrQT )[cosh(7TTQf) — 


y(CT,r;r^) := 


yR- 


-1 


= sinh(7rrg 


sin(7r<T)[cosh(7TTQf ) — cos 



Uqili := nali + kt'^tIz , s(<7, f ; t' ) — sinh(7TTg)[cosh(7rTof ) — cos(7rtr)] 1 

s(cr, f; Tq) = 7r / s(u, f; r^duli + 7ttq / s(a, u; Tg)dMl 2 . 
Jo Jo 

The differential equation for to solve in a G [— f , +1], f 6 [—1, +1], t > 0, is 



* f 4 +2* f2(f 2^ 2 +*^4^ 4 



The functional dependence of the physical measures, on the scaled variables and the param- 
eters, in q, is 



^(a,T,t;s;Re s ;Tg) , U(ct, f, t;s; Re s ; t' q ) , u)(a, f, f;s; i?e 5 ; r„) , P(ct, f, i; s; i?e 5 ; Tq) 
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The boundary and initial conditions for the streamfunction are 

*(flr,±l,t;5;ite J ,;7^)=TV(a;5;T^) , ±1, i?e s ; t' q ) = 

*(0,0,i;s;i?e s ;^) = , *(cr, f, 0; 5; i?e s ; t' ) = i>(a, f ;s; Re s ; t£) , 

where F(<7;s;to) = y(<J, ±1; Tq)s _1 = y'(<r; t^s -1 and ^(tr, f;s; i?e s ; Tq) denotes the initial 
flow field. 

Remark on the no-slip boundary condition: The motion of the cylinders is along the 
ii-direction. Therefore, the no-slip boundary condition is 

^ £ {a,±l,t,s;Re s ;T^) = . 

But, a = <t(i/';tq) and Tq £ = 0, therefore, 

[v{ti';To),f,t;s;Re s ;To] = ^a^y'Vi + *f f £ ■ 

For arbitrary values of f, y~ ^ 0. However, on the boundaries of the cylinders, it is 

dy'(p]T^) = djf(g,±l;r^) = 
dx(cr,±l;rg) dx(a, ±1; Tq) 

Therefore, the no-slip boundary condition is imposed by the requirement, that 

* f (CT,±l,t;s,.Re*;To) = . 

An asymptotic series solution for the streamfunction, valid for < Tq 2 <C 1 and 
< Tq Ren 1 is obtained by expressing the streamfunction and its boundary and initial 
conditions as 



f , t; s; Re s ; r ) = f , t; s)(r^e s )V^ j 

^(ct, f;s; Re,; t' q ) = \^{a, T;l)(r' Q Re s ) l T^' J 
J*(cr,±l,t;5) = TV(6r;7^;5)(5yo 
J*(0,0,t;5) = 



, ?* f (er,±l,t;s)=0 

, J*(er,f,0;5)=^(er,f;5) 



i5 denotes the Kronecker delta function. 

The solution for the streamfunction is obtained by substituting the assumed asymp- 
totic series in its differential equation and solving according to regular perturbation theory. 
The velocity and the vorticity are obtained similarly by the kinematic relations between the 
velocity or the vorticity and the streamfunction. 



*(a,f;i?e K ;r ) 



+ O{T' 2 Re R ) 



^3 



f)y' + 7 - 3f 3 + 2f)y%(^Re R )+ 



r^ + O(r 4 ) ; 



2' J " ' "7! 
(f 5 - 2f 3 + f)yi 2 + 0(4Re R ) 



U(a,f;s;i?e s ;r ) = tt" 1 [j* f (cr,f;5)^ 4i f;i)I 2 ] Ki?e 5 )V^' , 

2 V ■ 7 , 



U(ct,t;s; Re R ;T' ) = 



3 (f 2 - 1) y's- 1 + 3^ (7f 6 - 9f 2 + 2) y'y'-l-\r' Re R )+ 



+0(r^ 2 i? e 2 ? j] tt- 1 ^ 1 +0(7^)}li - 
+3^(f 7 - 3f 3 + 2f)[(^) 2 + y'y'- 2 ]l-^-\r' Re R ) + 0(rf Re R ) + 0{r' a 2 ) 



12 ; 
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ui(a,f;s;Re R ;To) = - 
+ O^ife^)! (^r^ 2 )- 1 + 0(1) 



For < Tq <C 1 and < Tg Rcr -c 1, when the asymptotic series solution is valid, 
the solution is independent of the initial condition and unique. 

For comparison, consider the case of irrotational fluid motion. For this case, the 
equation and the boundary conditions for the streamfunction iff (a, f ;s; t' ), are 

t^ 2 V 2 * = * f 2 + r^ 2 * 5 2 = , 

*(5,il;iy o ) = T7(ff;iy ) , *(0, 0; s; r£) = . 

An asymptotic series solution for the streamfunction, for irrotational flow, valid for < 
Tq 2 <C 1, is obtained by expressing the streamfunction and its boundary conditions as 

*(a,f;5;^)=^*(cr,f;i)^ 2j , j £ N , 
^(a,±l;B) = TV{a;8-y )5 J0 , 'tt(0,0;i)=0 . 

The solution for the streamfunction, for irrotational flow, is obtained by substituting the 
assumed asymptotic series in the Laplace equation for the streamfunction and solving ac- 
cording to regular perturbation theory to obtain 

f ;s; 4) = (-l) j+1 G(f; 2j)V^ rg 3 ' , j G No , where 
G(f; 2j) denote the polynomial functions of f , expressed recursively as 

G(f;0):=f ; G(f ; 2j - 1) := T J* + ~J ,VjeNi ; G(f; 2) := G(f; 1) ; 

G(f ; 2j) := G(f ; 2j - 1) - ^l^t iji , G No : fc £ [0, j - 2] , 
Vj 6 N : j e [2, +oo) . 
The velocity vector, for irrotational flow is then 

U(<7,f;s;^) = {-iy +1 G f {f;2j)V^- l T' 2j -% + (-l)'G(f ; 2j)V- a , i+ ^- 1 T ,2j l 2 . 

Remark on the unsteady fluid motion: The flow field, induced by the motion of 
the circular cylinders, is time dependent. The present formulation separates the effects of 
the unsteady fluid motion into two parts; a) fluid motion induced solely by and which is 
synchronised with the motion of the boundaries of the flow field and b) time-dependent 
fluid motion, which is independent of the motion of the boundaries. For irrotational flow, 
the second kind of fluid motion can be induced only by imposing time-dependent kinematic 
boundary conditions. For rotational flow, however, inertial effects, which are described by 
the di derivative, can also induce unsteady fluid motion of the second kind. 

For the present case of fluid motion, induced by the motion of the circular cylinders 
with constant velocity towards impact, the unsteady fluid motion is only of the first kind 
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at the asymptotic limit of approach of the cylinders. For irrotational flow, this is only 
because the kinematic boundary conditions are independent of time. For rotational flow, 
this is because of the steady boundary conditions and the fact that the c\ derivative does 
not appear in the differential equation for the leading order term of the asymptotic series for 
the streamfunction. The implication of this, is that the initial condition becomes redundant 
at the asymptotic limit considered and that the asymptotic series solution is independent of 
any initial flow field and thus, unique. 

5 An interpretation of the asymptotic series solution 

Consider the asymptotic series of the streamfunction for rotational and irrotational flow, 
respectively, 

^(a,f;Re R ;4)=i^(a,f)(4Re R ) i T^ j , fcfof;^) = j ^(a,f) T ^ j . 

The streamfunction for irrotational flow, is decomposed into infinite kinematic scales of 
volume flow rates of order 0[(U s R)tq 2: '), described by ^{a, f) . The streamfunction for 
rotational flow is also decomposed, into infinite kinematic scales of volume flow rates of 
order 0[{U s R)t^}. However, for this case, every kinematic scale is decomposed further, into 
infinite dynamic scales according to the dynamic scaling (T^Reu) 1 . The largest scale of the 
flow field of order 0(U s R) is described by ^(er, f). The lower order terms, jW(cr, f )(1 — Sijo), 
of the asymptotic series solution, are present and non-zero for all Tq, but, as Tq — > 0, 
they represent asymptotically vanishing physical flow scales of order O^sR^RcrYt^}. 
i^&(cr, f) are inter-dependent according to the differential equation 

^[^(a,f)](r^ efl )V^'=0, 

where 9 = (L> t ~V 2 - i?e- 1 V 4 ). 

According to the present interpretation, the initial condition ip(a, f ; Ren; Tq) = 
— {^(a, f)(roi?e_R)VQ 2 "', for the rotational flow field, in section HJ represents a multi-scale 
initial flow field. The magnitude of these scales, however, is adjusted according to the change 
of the Tq parameter, although the flow field at every scale can be arbitrary. 

At the limit Re R — > 0, the asymptotic series solution, for rotational fluid motion, 
becomes an asymptotic series solution for V 4 ^ = 0, as expected. 
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A Calculation of the basic differential operators in the 
s-space. 

The gradient V of a scalar /(s) is 

V/( S ) = / Si (s)I, . 
Si are independent variables; therefore, 

It = Im(i+1) x Im(i+2) = = Vs m ( j+ i) X Vs m ( i+2 ) . 

The divergence V- of a vector f (s) = /i(s)Ij is 

V • [/i(s)Ii] = V • {/i(s) • [Vs m(i+ i) x Vs„ l(l+2) ]} = 
= /t(s)V • [Vs m(l+1) x Vs m(i+2) ] + [Vs m(i+ i) x Vs m(j+2) ] • V/j(s) = 
= .A( s ){Vs m(j+2) • [V x Vs m(i+1) ] - Vs m(m) • [V x Vs m(i+2) ]}+ 
+ [Vs m(i+1) x Vs ro(i+2) ] • V/i(s) = I; • V/;(s) = / isi (s) . 

Therefore, 

V-f(s) = / i . < (s) . 
The local rotation of an arbitrary vector field V x f (s) is obtained as 

V x [Ms%] = V X [/i(s)VSi] = /i(s)V X Vs, + V/;( S ) X V Sl = 
= V/i(s) X Ij = fi 3m(i+2) (s)I m(i+ i) - fis m(i+1) (s)I ro (i+2) • 



Therefore, 



V X f(s) = [/m(t+2) Sm(i+i) (s) - /m(t+l) 5m(i+2) (s) 



I, • 



The Laplace operator V 2 of a scalar /(s) and the operator V 2 of a vector f (s) satisfies the 
identities 

V 2 /(s) = V • [V/(s)] , V 2 f (s) = V[V • f (s)] - V x [V x f (s)] . 

Therefore, 

V 2 /(s) = / s? (s) , V 2 f(s) = V 2 / 4 (s)I J . 

The biharmonic operator of a scalar /(s) and of a vector f (s), respectively, is 

V 4 /(s) = V 2 [V 2 /(s)] = f st (s) + / 5?si _» + / s?5i _ (s) 
V 4 f(s) = V 4 .A(s)I, . 

Therefore, in operational form, 



9 S .- Ii 



v 2 := a 



V 2 := J V 2 I, 



Vx 



l ( J +2) a 



Ii , 



9,4 + 9,2 2 + 9,2 2 
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